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Abstract. In this study we apply post-Newtonian (T-approximants) and 
resummed post-Newtonian (P-approximants) to the case of a test-particle in 
equatorial orbit around a Kerr black hole. We compare the two approximants 
by measuring their effectualness (i.e. larger overlaps with the exact signal), and 
faithfulness (i.e. smaller biases while measuring the parameters of the signal) with 
the exact (numerical) waveforms. We find that in the case of prograde orbits, T- 
approximant templates obtain an effectualness of ~ 0.99 for spins q < 0.75. For 
0.75 < q < 0.95, the effectualness drops to about 0.82. The P-approximants 
achieve effectualness of > 0.99 for all spins up to q = 0.95. The bias in the 
estimation of parameters is much lower in the case of P-approximants than T- 
approximants. We find that P-approximants are both effectual and faithful and 
should be more effective than T-approximants as a detection template family when 
q > 0. For q < both T- and P-approximants perform equally well so that either 
of them could be used as a detection template family. However, for parameter 
estimation, the P-approximant templates still outperforms the T-approximants. 
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1. Introduction 

Stellar mass compact binaries consisting of double neutron stars (NS), double black 
holes (BH) or a mixed binary consisting of a neutron star and a black hole, are 
the primary targets for a direct first detection of gravitational waves (GW) by 
interferometric detectors, LIGO 1 , VIRGO g], GEO600 0, and TAMA gj. Under 
radiation reaction the orbit of a binary slowly decays, emitting a signal whose 
amplitude and frequency increases with time and is termed a "chirp" signal. While it 
is believed that there is a greater population of NS-NS binaries El HI El OH , it is the 
BH-BH binaries that are the strongest candidates for detection since they can be seen 
from a greater volume, about two orders-of-magnitude greater than NS-NS binaries 

In order to detect such sources one employs the method of matched filtering 
The effectiveness of matched filtering depends on how well the phase evolution 
of the waveform is known. Even tiny instantaneous differences, as low as one part 
in 10 3 in the phase of the true signal that might be present in the detector output 
and the template that is used to dig it out could lead to a cumulative difference 
of several radians since one integrates over several hundreds to several thousands of 
cycles. In view of improving the signal-to-noise ratio for inspiral events there has been 
a world-wide effort in accurately computing the dynamics of a compact binary and the 
waveform it emits or to use phenomenologically defined detection template families 

There have been parallel efforts on using two different approximation schemes: 
On the one hand the post-Newtonian (PN) expansion of Einstein's equations has been 
used to treat the dynamics of two bodies of comparable masses with and without 
spin, in orbit around each other. This approximation is applicable when the velocities 
involved in the system are small but there is no restriction on the ratio of the masses 
[HO EH E| ED III EDI EH On the other hand, black hole perturbation theory has 
been used to compute the dynamics of a test particle in orbit around a spin-less or 
spinning black hole. Black hole perturbation theory does not make any assumptions 
on the velocity of the components, but is valid only in the limit when the mass of one 
of the bodies is much less than the other [221 EH1 EH ESI 1261 E3 • 

The post-Newtonian approximation is a perturbative method which expands the 
equations of motion, binding energy and GW flux as a power series in v/c, where v is a 
typical velocity in the system and c is the speed of light. At present, the PN expansion 
for the case of comparable-masses is known to order O (v 6 ) [20] and O (v 7 ) [21], for 
the energy and flux functions, respectively. 

As previously stated, black hole perturbation theory makes no assumptions about 
the orbital velocity of the components, but does restrict their masses. One assumes 
that a test particle of mass /x is in orbit about a central BH of mass M such that 
(i « M. Assuming this restriction is satisfied we have an analytical expression for the 
energy. However, no analytical expression has been worked out for the gravitational 
wave flux emitted by such a system. Using black hole perturbation theory, the 
most recent series approximation was calculated to O (v 11 ) by Tagoshi, Tanaka and 
Sasaki [23 ■ For a test particle in circular orbit about a Kerr black hole, the most 
recent progress is a series approximation to O (v 8 ) , by Tagoshi, Tanaka, Shibata and 
Sasaki [2211221 

Several authors [23 EDI ED E21 EH| have shown that the convergence of both post- 
Newtonian approximation and black hole perturbation theory is too slow to be useful 
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in constructing accurate templates. More recently, Damour, Iyer and Sathyaprakash 
(hereafter DIS) showed for the case of a test-mass in orbit about a Schwarzschild BH, 
that by using properly defined energy and flux functions that have better analytical 
properties, combined with Pade techniques, it was possible to take the existing series 
expansion and improve its convergence properties |33| . The new approximation in 
which Pade approximants of new energy and flux functions are used to derive improved 
templates is called P-approximant. While in general, more templates are needed for 
P-approximant templates to cover the same volume of parameter space |o 1| , the extra 
computational cost is preferred for the increased performance in P-approximants. 

In this paper we will extend the P-approximant technique to the case of a test 
particle orbiting a Kerr black hole. The reason for focusing on test-mass systems is 
that we can use the exact numerical fluxes jSSl from black hole perturbation theory 
with which to compare our results and thereby reliably demonstrate the usefulness of 
the technique. 

2. The Gravitational Waveform. 

In the stationary phase approximation the Fourier transform for positive frequencies 
reads 03 EH1 123 

h(f) ee ^ h(t) exp(27rz/i) dt = :^^L e *[W)-S] , (1) 

\ f 

where m = mi + m^, rj = mirr^/m 2 , C is a constant amplitude coefficient, d is the 
distance to the source, and, since h(t) is real, h(—f) = h*(f). Also, v — {itmf) 1 / 3 , / 
is the time-derivative of the instantaneous gravitational wave frequency evaluated at 
the stationary point given by, 

ttui hi'yv) 

where F is the gravitational wave flux function, E'(v) is the derivative of the orbital 
energy function with respect to v, i.e. E'(v) = dE/dv, and the evolution of the phase 
of the Fourier transform in the stationary phase approximation is given by solving 
following set of coupled I s * order differential equations 

dtp dt irm 2 E'(f) 



3. Gravitational Binding Energy and Flux Functions 

We can see from Equation © that the phase of the gravitational wave depends both 
on the energy and flux functions of the binary system. In the test-mass case we 
have an exact expression for the energy E(v), but only a series representation for 
the flux F(v). Numerically, for the test mass case, the flux has been computed 
exactly. Since our aim is to draw conclusions on how effective P-approximants are 
in the comparable mass case, wherein one has only a Taylor expansion of the flux, 
we construct P-approximants of the flux and compare it with the numerical results. 
Pade approximation can be thought of as an operator Pj^ that acts on a polynomial 

Sfc=o a k vk to define a rational function P$ = Y,k=o A k v k / 

(! + Efeii B k v k ^j such 

that the N + M + 1 coefficients in the rational polynomial on the right hand side is 
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the same as the n+l Taylor coefficients on the left hand side. By setting N = M + e 
with e = 0, 1, we can define two types of Pade approximants: These are the super- 
diagonal, Pm +<l , and sub-diagonal, Pj$ +e , approximants. Normally, the sub-diagonal 
approximants are preferred over super-diagonal approximants. This is because when 
M = N + e the rational function can be re-expanded as a continued fraction which 
has the property that as we go to each new order of the power series only one new 
coefficient needs to be calculated. Conversely, with the super-diagonal approximants, 
we would have to re-calculate all the A's and B's in the above equation as we go 
to higher orders in the Taylor expansion. This means that the sub-diagonal Pade 
approximants are more stable and if we see a trend of convergence in the coefficients 
the addition of a term is not likely to spoil this convergence. 



3.1. The Orbital Energy. 

In the case of both Schwarzschild and Kerr black holes we have an exact expression 
for the orbital energy of a test particle in a circular orbit around the parent black 
hole. For a black hole of mass M the energy E in terms of the dimensionless 
magnitude of velocity v = <JM/r, r being is the radial coordinate in the Boyer- 
Lindquist coordinates, takes the form 001 

y/l - 3v 2 + 2qv 3 

where q is a dimensionless spin parameter given in terms of the spin angular 
momentum J of the black hole by q = J /M 2 = a/M, with a spin angular momentum 
per unit mass in the Kerr metric. 



3. 2. Post-Newtonian flux function 

For a test-particle in a circular equatorial orbit, the post-Newtonian expansion of the 
flux function has been calculated up to O (w 11 ) in the case of a Schwarzschild BH |27j. 
and to 0(v s ) in the case of a Kerr BH |281 129| . The general form of the flux function 
in both these cases is given by the expression 



F Tn (x; q)=F N {x) 



]T a k (q)x k + ln(x) h(q)x k + O (x n+1 ) 

,fe=0 fc=6 



,(5) 



where F^(x) — ^r-r) 2 x 10 is the dominant Newtonian flux function. Here, x is the 
magnitude of the invariant velocity parameter observed at infinity which is related to 

1 /3 

the angular frequency f2 by x = (Mil) . The relation between the parameter x and 
the local linear speed v in Boyer-Lindquist coordinates is given by 

x(v, q) = v [1 - qv 3 + q 2 v 6 ] 1/3 , (6) 
which reduces, in the Schwarzschild limit, to x = v. 



3.3. P-approximant of the flux function 



In order to prepare the series representation of the flux for creating the Pade 
approximation, it is convenient if we factor out the logarithmic terms. We can then 
write Equation (JSJ as 



F Tn (x) =F n (x) 



1 + ln 



x 

•Else 



fe=6 



J Lfe=0 



(7) 
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where the new coefficients Cfe and Ik are functions of the old coefficients au and bk . As 
in Reference [23 the log-terms have been "normalized" using the value of the velocity 
parameter at the LSO; this helps in reducing the importance of the log-terms. xi so is 
found by substituting vi so = \J M/ri so into Equation JfJJ, where |4(Jj 



rL(<l)= M 3 + Z 2 {q) T >/[3 " Zi(q)] [3 + Zl (q) + 2z 2 (q)} 



where 

Zl (q) = 1 + (1 - q 2 ) 



(l + q)' 3 +{l-q)- 



z 2(q) = \/3<7 2 + z 



(8) 



(9) 



and the + (— ) sign corresponds to prograde (retrograde) orbits. We create the factored 
flux function, /t„(x) by the operation fx n = (1 — x/x po i e ) Fx n , where x po i e is obtained 
in the same way as xi so but by using llT] 



pole 



(q)=2M 



1 



cos 



- cos 1 (tq) 



(10) 



Factoring out the pole also helps to alleviate the problem arising from the absence 
of the linear term in the PN expansion of the flux. (Note that oi = in both the 
Schwarzschild and Kerr cases.) If we write the expression in full we obtain 



f Tn (x)=F N (x) 



1 + ln 



x 

%lso 



E l ^ 

k=6 



E 

,fe=0 



(ii) 



where f = c g and f k =c k ~- c k _Jx pole , k=l,...,n, 

We can now construct a new flux function by constructing the polynomial 
expansion of the inverse of the flux function and construct the Pade approximant 
of the resulting polynomial. We call the approximant constructed this way as Inverse- 
cm I-Pade approximant because it is obtained from the Taylor expansion of the inverse 
of the flux function f(y) in Equation The Inverse Pade approximant of flux is 

defined by 



f IPn (x)=F^(x) 



1 - In 



x 

•Else 



k=6 



E<4 

,fe=0 



(12) 



where the coefficients in the Taylor expansion are defined by X)fc=o d k % k = 
{^2k =Q f k x k ) 1 ■ We call the flux function constructed in this manner P-approximant. 
Thus, we define the Inverse P-approximant as, 



Fip„(x) = 



1 - 



%pole 



(13) 



4. Effectualness and Faithfulness of T- and P-Approximants 

We shall address the performance of the approximants in extracting the exact 
waveform in two ways: The effectualness of the templates measured in terms of the 
maximum overlap they can achieve with the exact waveform when the parameters of 
the approximant are varied in order to achieve a good match. The faithfulness of the 
approximant templates measured in terms of the systematic errors in the estimation 
of parameters while detecting exact waveforms. 
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Figure 1. The maximized prograde overlaps for T-approximant (left) and P- 
approximant (right) templates at the x s approximation. Each system consists of 
a 1.4 Mq NS inspiralling into a central BH of mass ranging from 10-50 Mq. 



4-1- Overlaps and fitting factor 

We define the scalar product of two waveforms h and g by 

df 



(h\g) 



s h (f) 



h(f)f(f) + h*(f)~g(f) 



(14) 



where the * denotes complex conjugate and h(f), g(f) are the Fourier transforms of 
h(t), g(t). For initial LIGO, the one-sided noise power spectral density (PSD) from 
the design study is given by [HHI 



S h (f) = 9x 1CT 46 [0.52 + 0.16a;- 4 - 52 + 0.32x 2 ] Hz~ 



(15) 



where x = f/fk, and fk — 150 Hz is the "knee-frequency" of the detector. We take 
the PSD to be infinite below the lower frequency cutoff of f\ ow = 40 Hz. For two 
normalized waveforms, or signal vectors, the scalar product returns the cosine of the 
angle between them and is normally referred to as the overlap. For detection of signals 
what is more important is the fitting factor FF : As each template is a function of 
the parameters A* 1 , the fitting factor is defined as the maximum overlap obtained by 
varying the parameters of the template (or the approximate waveform) relative to the 
exact waveform: 



FF 



max O (A M ) 



(16) 



If two waveforms are a perfect match then their overlap is unity. In the case 
of circular equatorial orbits of a test mass around a central black hole, the wave is 
parameterized by the set A M = (to, m i l)- Maximization over t is achieved 
by simply computing the correlation of the template with the data in the frequency 
domain followed by the inverse Fourier transform. It was pointed out |42U43| that the 
maximum of the overlap of the data with a template over $o can be computed using 
just two templates - an in-phase and a quadrature-phase template. 

Using a zero-phase un-normalized waveform, ho = h ($o = 0), we generate two 
orthonormalized waveforms according to Hq = ho/\ho\ and H„/2 — iHo, which is 
explicitly orthogonal to the in-phase template. The (square of the) maximum of the 
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" (V, q) x (v, q) 

Figure 2. The 4-PN T-approximant (left) and P-approximant (right) analytical 
flux function for spins of q = 0.5, 0.75 and 0.95 against the numerical fluxes for 
the same spin values. 

overlap over <£>o is given by the sum-of-squares of the overlap of the signal with the 
in-phase and quadrature-phase templates: 

max (0) = J{H \h*) 2 + (H« /2 \h*) 2 , (17) 

where h x denotes the "exact" waveform. Once this is done, we use a maximization 
routine to find the optimal values of the parameters (m, rj, q). 

In this study, as we are working in the test-mass approximation, we assume that 
our system is composed of objects with a small mass ratio. For concreteness we assume 
that the system comprises a 1AM@ NS inspiralling into a central BH of varying spin 
magnitude and mass. Beginning with a central BH of IOMq, we work our way upward 
to a 50A/q BH. We will also look at the limiting case of a 10-10M Q equal-mass system. 
Strictly speaking, the formulas for energy and flux functions used in this study are 
not applicable to the comparable mass case since we have neglected the finite mass 
correction terms in these quantities. However, the results of such a study should give 
us an indication of how strong are the relativistic corrections, as opposed to finite- 
mass corrections, in the case of comparable masses. In all cases we vary the spin 
magnitude of the central BH from q — —0.95 to 0.95. Since the main focus in this 
study is test-mass approximation we shall be interested in errors in estimating the 
chirp- mass, A4 — m?/ 3 / 5 , in addition to the spin magnitude of the central BH. In 
view of economy we shall only present the results for the highest PN order available, 
namely 0(x 8 ) order. In all cases our fiducial exact signal h x will be that obtained 
by using the exact expression for the energy in Equation and the exact numerical 
fluxes generated using black hole perturbation theory |35j . and the template will be 
the approximate waveform constructed using the exact expression for the energy, as 
before, and an approximate expression for the flux. 

4-2. Prograde Orbits - Effectualness 

For T-approximants - Figure ^ left panel - maximizing over all parameters gives 
fitting factors of FF > 0.98 at all mass ranges for the test-mass systems up to a spin 
of q = 0.75. Between 0.75 and q — 0.95, the T-approximant templates begin to perform 
badly and the fitting factors drop to FF ~ 0.82. For the equal-mass case - Figure0- 
the templates once again achieve fitting factors of FF > 0.99 up to a spin of q = 0.75, 
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Figure 3. The percentage error in the estimation of the chirp-mass, M, for T- 
approximant (left) and P-approximant (right) templates at the x 8 approximation. 

but fall off at higher spin magnitudes achieving a fitting factor of ~ 0.98 at q = 0.95. 
We should point out that these results do not properly convey just how bad the 4-PN 
T-approximant template actually performs. In the left hand panel of Figure [2 we 
plot the PN approximation for the flux function against the numerical fluxes at spins 
of q = 0.5,0.75 and 0.95. These have corresponding values of xi so w 0.48,0.54 and 
0.66. We can see that the flux function at q — 0.75 and 0.95 become negative long 
before the LSO is reached. This means that as we go to higher and higher spins, we 
can model less and less of the waveform. For q < 0.6 we can model the waveform up 
to the LSO. However, for q > 0.6 this we have to stop the waveform generation at a 
cutoff velocity of x cut « 0.5. Therefore, the fitting factors beyond q = 0.75 completely 
overestimate the performance of the T-approximant template. 

If we now focus on the right hand panel of Figure [3 we see the true power of the 
P-approximant templates. The 4-PN template suffers none of the divergences that 
effect the T-approximants. We therefore generate all templates up to the LSO or 2 
kHz, whichever is reached first. We can see from the right panel of Figure ^ that 
the P-approximant templates achieve fitting factors of > 0.99 at all spin values. For 
the equal-mass case - Figure [5]- the P-approximant templates achieve fitting factors 
of FF > 0.995 at all mass and spin levels. This demonstrates that in the case of 
prograde orbits, the P-approximant templates are clearly more robust, even at high 
spin magnitudes of q = 0.95. 

4-3. Prograde Orbits - Faithfulness 

In Figure |3| we have plotted the percentage bias in the estimation of the chirp-mass 
for both T- (left panel) and P-approximant (right panel) templates. From Figure [21 
left panel, it is clear that for T-approximants the bias in M. varies between - 5%. 
Comparing the left and right panels of Figure |3 we find that the P-approximant 
templates (right panels) are more faithful with the percentage bias being less than 1% 
in general. 

In Figure0]we present the percentage bias in the estimation of the spin magnitude 
of the central BH. We find that for both approximants we have to endure a large bias 
in q. Once again the P-approximant templates are more faithful. For T-approximants 
the bias varies between 5 and 90%, while for P-approximants it lies between 2 and 10%. 
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Figure 4. The percentage error in the estimation of the spin parameter, q, for T- 
approximant (left) and P-approximant (right) templates at the x 8 approximation. 

From the middle and bottom Sections of Figure |S] it is clear that even in the equal 
mass case the bias in the estimation of both parameters is greater for T-approximant 
templates than P-approximants. It is therefore clear that in the case of prograde 
orbits, P-approximant templates must be used in any detection strategy. 

4-4- Retrograde Orbits - Effectualness 

In the case of retrograde orbits both templates achieve fitting factors of FF > 0.99 
regardless of the spin magnitude and the chirp mass. There is no surprise here: The 
retrograde waveforms are still well within the adiabatic regime and are, therefore, 
modelled well by both templates. From a purely detection point of view, unlike the 
prograde case, there is no obvious benefit from employing P-approximant templates. 
For the equal mass system too there is not much difference in the fitting factors of the 
two families of templates with the exact signal. 

4-5. Retrograde Orbits - Faithfulness 

The benefit of using P-approximant templates for retrograde motion is only observed 
when we consider parameter extraction. Referring to the bottom portions of the two 
panels in Figure |3 we note that the T-approximants perform well at all spins for 
3.0 < A4/Mq < 4.5, with a bias of less than 1% in the estimation of Ai. Beyond 
this, there is in general a bias of > 2%. As we approach the extreme retrograde case 
the bias rises to as much as 55%. The P-approximants perform in a similar manner. 
The bias in the region 3 < M/Mq < 4.5 is again in general less than 1%. The error 
does again increase as we head towards the extreme test-mass range, but in this case 
it reaches a maximum value of 8% as opposed to the 55% seen in the case of the 
T-approximants. 

The error in estimating q reaches a maximum of 80% for T-approximants and 25% 
for P-approximants. This is consistent with the results presented in Reference |28l 129] 
where it was shown that the PN approximation waveforms perform worse in the 
retrograde case. Now, referring once more to the equal mass case, the P-approximant 
templates outperform their T-approximant counterparts. Even with this, we must 
again conclude that while on the surface there is no clear case for using P-approximant 
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Figure 5. The fitting factors for a 10-10 Mq binary - without the finite mass 
correction terms - for T and P-approximant templates (top) . The percentage error 
in the estimation of the chirp-mass for T and P-approximant templates (middle). 
The percentage error in the estimation of the spin of the central black hole for T 
and P-approximant templates (bottom). 



templates in searching for retrograde motion systems, they should be used because of 
the lower bias incurred in the estimation of parameters. 

5. Conclusions 

We have applied P-approximant templates to the case of a NS orbiting Kerr BHs of 
varying mass and spin. Using a signal waveform constructed from the exact expression 
for the orbital energy and numerical fluxes from black hole perturbation theory, we 
were able to compare the performance of T and P-approximant templates. In the case 
of retrograde, Schwarzschild and prograde orbits, not only did the P-approximants 
gave better and more reliable fitting factors, they also gave smaller biases in the 
estimation of parameters. We also saw the true power of the P-approximant templates 
in that we were able to generate templates right up to the LSO . This is something 
that was not possible with the T-approximant templates due to the approximation 
for the flux function becoming negative before the LSO is reached. While not being 
completely correct due to the fact that we omitted the finite-mass correction terms, 
we also saw that the P-approximant templates gave the best performance when in the 
equal-mass case. 

It is clear that for the type of systems examined in this paper, namely equatorial 
test-mass circular orbits in Kerr, P-approximant templates are to be preferred over 
their PN counterparts in both detection and measurement for prograde systems. For 
retrograde systems, while the T-approximants can be used for detection purposes, the 
P-approximants are superior when it comes to parameter estimation. 
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